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Abstract 

We estimate the sensitivity to nano-Hertz gravitational waves of pulsar timing experiments in 
which two highly-stable millisecond pulsars are tracked simultaneously with two neighboring radio 
telescopes that are referenced to the same time-keeping subsystem (i.e. "the clock"). By taking 
the difference of the two time-of-arrival residual data streams we can exactly cancel the clock noise 
in the combined data set, thereby enhancing the sensitivity to gravitational waves. We estimate 
that, in the band (10 -9 — 10 -8 ) Hz, this "interferometric" pulsar timing technique can potentially 
improve the sensitivity to gravitational radiation by almost two orders of magnitude over that of 
single-telescopes. Interferometric pulsar timing experiments could be performed with neighboring 
pairs of antennas of the forthcoming large arraying projects. 
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The direct detection of gravitational waves is one of the most challenging experimental 
efforts in physics today. A successful observation will not only represent a great triumph in 
experimental physics, but will also provide a new observational tool for obtaining a better 
and deeper understanding about their sources, as well as a unique test of the various proposed 
relativistic metric theories of gravity pp. 

Pulsar timing experiments aimed at the detection of gravitational radiation have been 
performed for decades now. The basic principle underlining the pulsar timing technique is 
the same as that of other gravitational wave detector designs [2HIj: to monitor the frequency 
variations of a coherent electromagnetic signal exchanged by two or more "point particles" 
separated in space. As a pulsar continuously emits a series of radio pulses that are received 
at Earth, a gravitational wave passing across the pulsar-Earth link introduces fluctuations in 
the time-of-arrival (TOA) of the received electromagnetic pulses. By comparing the pulses 
TOAs against those predicted by a model, it is in principle possible to detect the effects 
induced by any time- variable gravitational fields present, such as the transverse-traceless 
metric curvature of a passing plane gravitational wave train [8]. The frequency band in 
which the pulsar timing technique is most sensitive to ranges from about 1CT 9 to 10 -6 
Hz, with the lower-limit essentially determined by the overall duration of the experiment 
(1CT 9 Hz ~ 1/30 years), and the upper limit identified by the signal-to-noise ratio of the 
received radio pulses. To attempt observations of gravitational waves in this way, it is thus 
necessary to control, monitor and minimize the effects of other sources of timing fluctuations, 
and, in the data analysis, to use optimal algorithms based on the different characteristics 
of the pulsar timing response to gravitational waves (the signal) and to the other sources of 
timing fluctuations (the noise). 

A quantitative analysis of the noise sources affecting millisecond pulsar timing searches 
for gravitational radiation has been highlighted in a recent paper by Jenet et al. [9]. There 
it has been shown that, in the region (1CT 9 — 10 -8 ) Hz of the accessible frequency band, the 
main noises are due to: 

(i) finiteness of the signal-to-noise ratio in the raw observations; 

(ii) uncertainties in solar system ephemerides, which are used to correct arrival times at 
the Earth to the barycenter of the solar system; 

(iii) variation of the index of refraction in the interstellar and interplanetary plasma; 
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(iv) intrinsic rotational stability of the pulsar, and 

(v) instability of the local clock against which pulsars are timed, and noise in time transfer 
if the clock is not located at the observatory site. 

Although the timing fluctuations induced by most of these noises can be in principle ei- 
ther reduced or calibrated out, the fundamental noise-limiting sensitivity of pulsar timing 
experiments is imposed by the timing-fluctuations inherent to the reference clocks that con- 
trol the TOA measurements. The noise analysis in [9| indicated that clocks such as the 
Linear-Ion- Trapped-Standard (LITS) (presently in-the-field state-of-the art atomic clocks 
with long-time- scale timing stability) would result in a sinusoidal strain sensitivity of pulsar 
timing searches for gravitational waves to a level of about 1CT 15 after coherently integrating 
the data for a period of 10 years. Since the characteristic wave amplitude associated with 
a super-massive black-hole-binaries background is predicted to be of comparable magnitude 
at the frequency 3 x 1CT 9 Hz pm [TT], it is clear that a single telescope will not be able to 
unambiguously detect such a gravitational wave signal. 

A method for statistically enhancing the signal-to-noise ratio of pulsar timing experiments 
to backgrounds of gravitational waves was first proposed by Hellings and Downs [12J and 
further analyzed and improved by Jenet et al. [13J. This technique relies on cross-correlating 
pairs of TOA residuals data taken by an array of radio telescopes observing highly-stable 
millisecond pulsars. Since the gravitational wave background is common to all the timing 
residuals while many noises affecting them are in principle uncorrelated, the cross-correlation 
technique enhances the strength of the gravitational wave signal over that of the noise. 

The experimental configuration we propose in this letter is complementary to the cross- 
correlation method, in that it provides a way for exactly canceling the clock noise from pairs 
of TOAs residuals. This is possible by relying on TOAs residuals measured simultaneously by 
two neighboring radio telescopes tracking, with the same clock, two highly-stable millisecond 
pulsars located in two different parts of the sky. Since the two TOAs responses to the 
gravitational wave signal will be different while the clock timing noise transfer function in 
the two TOA residuals will be the same, by taking the difference of the two time series 
the clock noise is exactly canceled from the resulting synthesized data while sensitivity to a 
gravitational wave signal is preserved. 

Let us denote with R\(t), R2(t) the two time-of- arrival residuals measured at time t at 
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the two telescopes. Their time derivatives can be written in the following way 

d ^ = (1 - k • [*i(t - u{\ + k • n,)) - %(t)} 

+C(t) + Ni(t) i = l,2, (1) 

where we have denoted with C(t) the random process associated with frequency fluctuations 
due to the clock (common to the two timing residual measurements), with iVj(t) , i = 1, 2 the 
two random processes associated with frequency fluctuations due to all other noise sources 
affecting the timing residuals, and the function ^i(t) contains the contribution of a spin-2 
gravitational wave transverse-traceless strain tensor h in the following way 

_ , . (hi ■ h ■ hi) 
2[l-(k-n,)2] 

In Eqs.([lJ [2J AjLj , z = 1,2 is the vector oriented from telescope i to pulsar i, k is the 
unit vector associated with the direction of propagation of the wave, and the symbol • is the 
operation of scalar product (see Fig. [I]). If we denote with AR the difference of the time- 
derivative of the two timing residuals given in Eq. (fl]), we find the resulting interferometric 
pulsar timing response 

AR{t) = R^t) - R 2 {t) = Nt(t) - N 2 (t) 

+ [(l-k-n 2 ) ^-(l-k-nO tf x (*)] 
+ [(l-k-h^xit-Lxil + k-hx)) 

- (l-k-n 2 )^ 2 (t-L 2 (l + k-n 2 ))] . (3) 

Note how the two-pulse time-structure of the time-derivative of the two TOA residuals 
given in Eq. ([!]) becomes a three-pulse response in the interferometric pulsar timing expres- 
sion given by Eq. ([3]); the three pulses correspond to the times a gravitational wave interacts 
with the Earth (which are the same in the two TOAs residuals) and the two millisecond 
pulsars. 

To assess the effectiveness of the interferometric pulsar-timing technique we have esti- 
mated its sensitivity by relying on the noise-model discussed in |9]. In particular, we have 
assumed that: (i) multiple-frequency measurements can be implemented in order to ade- 
quately calibrate timing fluctuations due to the intergalactic and interplanetary plasma; (ii) 
more accurate solar system ephemerides will be available by the time these experiments 



will be performed; and (iii) the tracked millisecond pulsars have frequency stabilities better 
than those of the operational ground clocks. A recent stability analysis of presently known 
millisecond pulsars [T7] has shown that there might exist some with frequency stabilities 
superior to those displayed by the most stable operational clocks in the (1CT 9 — 1CT 8 ) Hz 
frequency band. It should be emphasized, however, that the interferometric pulsar timing 
experiments discussed in this letter can also be used to isolate and study the intrinsic pulsars 
spin-noises with a much higher precision than what currently achievable by single-antenna 
experiments. 

The sensitivity of a detector of gravitational waves has been traditionally taken to be equal 
to (on average over the sky and polarization states) the strength of a sinusoidal gravitational 
wave required to achieve a given signal-to-noise ratio (SNR) over a specified integration time, 
as a function of Fourier frequency. Here we will assume a SNR = 1 over an integration time 
of 10 years. In order to compute the interferometric pulsar timing sensitivity we will consider 
two millisecond pulsars of comparable intrinsic stability and at a distance of 1 kpc from the 
Earth. We will also assume the following expression for the power spectrum of the noises 
Ni , % = 1, 2 introduced in Eq. Q 



which corresponds to a white-timing noise of 100 nsec. [9] in a Fourier band +/- 0.5 cy- 
cles/day (i.e. one sample per day). The 100 nsec. level is the current timing goal of leading 
timing array experiments as three pulsars are being timed to this level [17]. 

The gravitational wave sensitivity of our proposed interferometric pulsar timing technique 
is then defined as the wave amplitude required to achieve a signal-to-noise ratio of 1 in a 
ten years integration time: \J2Sn{J) 5/(root- mean-squared gravitational wave response). 
Note the factor 2 multiplying the spectrum Sn follows from treating the noises iVj , i = 
1,2 as independent; the bandwidth, B, was taken to be equal to one cycle/10 years (i.e. 
3.17 x 10~ 9 Hz). In order to calculate the root-mean-squared gravitational wave response 
we have assumed waves to be elliptically polarized and monochromatic, with their wave 
functions, (h + , h x ), written in terms of a nominal wave amplitude, H, and the two Poincare 
parameters, ($,T), in the following way [18] 



S N (f) = 3.4 x 10~ 8 f Hz 



-l 



(4) 



h x (t) = H cos(r) sin(u;t) . 



h+(t) = H sin(r) sm(ut + $) , 



(5) 
(6) 
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FIG. 1. Two millisecond pulsars are tracked simultaneously by two radio telescopes sited near 
each other. The unit vector k defines the direction of propagation of the gravitational wave signal 
and (9,(j)) are the usual polar angles associated with it. The two unit vectors, (111,112), define the 
directions in the sky to the two pulsars, and the angle between them has been denoted with 2Q. 
The bisector of this angle coincides with the x-axis of an orthogonal coordinate system in which 
the y-axis lies in the plane defined by the two unit-vectors (111,1x2), and the z-axis is chosen to 
form a right-handed orthogonal triplet (x, y, z). 

We averaged over source direction and polarization states by assuming uniform distribution 
of the sources over the celestial and polarization spheres respectively. The averaging was 
done via Monte Carlo integration with 10000 source position/polarization state pairs per 
Fourier frequency bin and 5 x 10 5 Fourier bins across the (10 -9 — 10~ 6 ) Hz band. In Figure 
(J2J) we plot the sensitivity of the interferometric pulsar timing technique as a function of the 
Fourier frequency and for various values of the angle 2Q enclosed between the directions to 
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FIG. 2. Sensitivity of the interferometric pulsar timing technique as a function of the Fourier 
frequency and for various values of the pulsars angular separation in the sky, 20, . The pulsars 
have been assumed to be at equal distance from Earth (1 kpc). Note the mild dependence of 
the sensitivity over the angle 20 in the interval of values considered. The sensitivity of a single- 
telescope experiment (blue curve), as well as the current estimates of the upper and lower bounds 
on the characteristic amplitude radiated by an ensemble of super-massive black-hole binaries (two 
black dashed-lines) are included. See text for details. 

the two pulsars. For sake of simplicity we have assumed the two pulsars to be at an equal 
distance of 1 kpc from Earth. We have independently analyzed configurations with unequal 
distances and found no significant differences to the sensitivity curves shown in Fig.(|2]). 
Another property of the sensitivity curves is their mild dependence over the angle 2Q; a 
configuration with two millisecond pulsars that are separated by only 15° in the sky shows a 
sensitivity that is only 28 percent worst than that with 2Q = 90°. However, for values of 2Q 
smaller than 15°, the sensitivity rapidly worsen and goes to zero at 2Q = 0. This is of course 
consequence of having assumed, for this calculation, the two pulsars to be at equal distance 
from Earth. In the case when L\ ^ L 2 instead, the sensitivity to gravitational radiation 
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does not go to zero at 2f2 = 0, as it can be easily understood by simple inspection of Eq. 
(j3|. Since longer joint tracking times can be achieved with pairs of millisecond pulsars that 
have smaller angular separation, our analysis indicates that these configurations (with, let's 
say 2Q ~ 20°) might be preferable to those with larger angle 2Q. 

For comparison, in Fig.(|2]) we have included the sensitivity of a single-telescope experi- 
ment that relies on a LITS clock [9]. Although at frequencies larger than about 10~ 8 Hz this 
is about a factor of 2 better than that of the interferometric technique (in this part of the 
band the noises affecting the interferometric combination are twice as many), at frequencies 
lower than 10~ 8 Hz the advantage of the interferometric technique becomes evident. In this 
region of the band the clock noise dominates all other noise sources in a single-telescope 
timing residuals data, and our technique exactly cancels it. 

In order to put in perspective the sensitivity enhancement brought by the interferomet- 
ric pulsar timing technique over single-telescope experiments, we have included the current 
estimated range of strengths of an astrophysical gravitational wave background from inco- 
herently radiating super-massive black hole binaries [TOl ITT] . At 3 x 10~ 9 Hz we find that 
the most conservative prediction of the strength of such a background would result into an 
SNR ~ 5, while a more optimistic estimate of its characteristic strain amplitude would make 
it detectable with a SNR ~ 50. 
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